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Background: Whether atomic nuclei can possess triaxial shapes at their ground states is still a subject of 
ongoing debate. According to theory, good prospects for low-spin triaxiality are in the neutron-rich Mo-Ru region. 
Recently, transition quadrupole moments in rotational bands of even-mass neutron-rich isotopes of molybdenum 
and ruthenium nuclei have been measured. The new data have provided a challenge for theoretical descriptions 
invoking stable triaxial deformations. 

Purpose: To understand experimental data on rotational bands in the neutron-rich Mo-Ru region, we carried out 
theoretical analysis of moments of inertia, shapes, and transition quadrupole moments of neutron-rich even-even 
nuclei around ^^°Ru using self-consistent mean-held and shell model techniques. 

Methods: To describe yrast structures in Mo and Ru isotopes, we use nuclear Density Functional Theory (DFT) 
with the optimized energy density functional UNEDFO. We also apply Triaxial Projected Shell Model (TPSM) 
to describe yrast and positive-parity, near-yrast band structures. 

Results: Our self-consistent DFT calculations predict triaxial ground-state deformations in ^°®’^°®Mo and 
i°® iiO’ii^Ru and reproduce the observed low-frequency behavior of moments of inertia. As the rotational fre¬ 
quency increases, a negative -7 structure, associated with the aligned i'(hii/ 2 )^ pair, becomes energetically favored. 
The computed transition quadrupole moments vary with angular momentum, which rehects deformation changes 
with rotation; those variations are consistent with experiment. The TPSM calculations explain the observed band 
structures assuming stable triaxial shapes. 

Conclusions: The structure of neutron-rich even-even nuclei around ^^°Ru is consistent with triaxial shape de¬ 
formations. Our DFT and TPSM frameworks provide a consistent and complementary description of experimental 
data. 

PACS numbers: 21.60.Jz, 21.60.Cs, 21.10.Re, 21.10.Ky, 27.60.+j 


I. INTRODUCTION 

Neutron-rich molybdenum and ruthenium isotopes are 
known to exhibit shape changes and shape-coexistence 
phenomena [IH6]. With increasing neutron number, tri¬ 
axial deformation is expected to appear in their ground 
states due to the occupation of lvhii/2 and l7rgg/2 in¬ 
truder orbitals [?]• 

Experimentally, the clearest signature of triaxial 
shapes comes from the 7-ray spectroscopy of rotating nu¬ 
clei BM- Following Ref. [ 10 ], which reported evidence 
for rotational-like behavior in the very neutron-rich even- 
even Zr-Pd region, numerous experiments were devoted 
to investigations of shape transitions and rotational prop¬ 
erties in this region. The first systematic high-spin study 
of collective band structures was undertaken in Ref. 
which reported deformed conhgurations in io3.io4.io72j. 
and 

The early work [ 5 ] on ^^"‘Ru was indicative of a tran¬ 
sition from spherical to triaxial shapes. The collective 
nature of neutron-rich nuclei and triaxiality in ^^*^Ru and 
^^^Ru was confirmed in Ref. (T^] by reporting a steady de¬ 
crease of the 7-band bandhead energy. More evidence of 


collective triaxial behavior of Ru isotopes came from the 
spectroscopy of fission fragments [TSHTS] . The decrease 
of transition quadrupole moments at high spin showed 
that the triaxial deformation in neutron-rich Mo isotopes 
could be spin-dependent m- Another piece of experi¬ 
mental information on triaxial deformations came from 
the measurement of the quasi-7 band in ^^°Mo ^01121] . 
A more detailed information on quadrupole collectivity 
was obtained by Coulomb excitation studies, which suc¬ 
ceeded in determining unique sets of E 2 and Ml matrix 
elements in ^°"^Ru | 3 ] and ^^°Mo [ 3 | and extracting triax¬ 
ial deformations using the collective quadrupole invariant 
approach. Rotational bands in ^°®Mo and i08,iio,ii2j^^ 
were investigated in Ref. [H [ 22 ] , which reported chiral 
doublets associated with triaxial nuclear rotation. In re¬ 
cent papers H 113], transition quadrupole moments of 
rotational bands in neutron-rich, even mass ^°^“^°®Mo 
and ^°®“^^^Ru nuclei were measured for the spin range 
of 8-16 S, suggesting 7-softness effects or even triaxiality 
in these nuclei. 

Theoretically, triaxial ground states in this region have 
been investigated with different models. In Ref. | 7 ], based 
on macroscopic-microscopic approach, triaxial ground- 
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state (g.s.) minima were found in the neutron-rich Mo 
isotopes with N=62-66 and also in Ru isotopes. In a sys¬ 
tematic survey of Refs. [211[5S], the largest shell effects 
due to triaxial deformations were found around ^"^^Ru. 
The interacting boson model analysis of Ref. [3^ did not 
find any candidates for stable triaxiality in this region. 
In Ref. [ 37 ], potential energy surface (PES) calculations 
for Ru isotopes were carried out with Hartree-Fock (HE) 
and interacting boson models, and shallow triaxial min¬ 
ima were found for N = 64 — 70 (see also Ref. [28]). 
In the self-consistent Hartree-Fock Bogoliubov study of 
Ref. [ 33 ] with Gogny DIS interaction, triaxial deforma¬ 
tions were predicted for even-even isotopes 
and 

Collective rotation has been shown to enhance triaxial 
minima in even-even Mo and Ru isotopes [HEOIISI]. For 
that reason, those nuclei are candidates for the presence 
of novel collective modes associated with triaxial rota¬ 
tion, such as wobbling motion and chiral bands mm- 
The angular momentum alignment pattern in the lowest 
bands of Mo and Ru isotopes was also explored within the 
projected shell model approach [331 [33] • The axial study 
[ 33 ] provided a reasonable description of yrast spectra 
and electromagnetic properties of ^°°“^^^Ru. This work 
was extended in Ref. [ 33 ], which also contains TPSM 
analysis of ^^°Mo and ^^^Ru. In the case of ^^"^Ru, stable 
7 deformation turned out to be crucial for reproducing 
the data. 

In an attempt to explain the recent data on transition 
quadrupole moments in Mo and Ru nuclei, cranked rel¬ 
ativistic mean-field calculations |23j predicted axial pro¬ 
late and oblate ground states in those nuclei. However, 
the angular momentum dependence of resulting transi¬ 
tion quadrupole moments was not consistent with obser¬ 
vations. As concluded in Ref. [ 33 ]: Attempts to describe 
the observations in mean-field based models, specifically 
cranked relativistic Hartree-Bogoliubov theory, illustrate 
the challenge theory faces and the difficulty to infer in¬ 
formation on '^-softness and triaxiality from the data. To 
shed some light on this puzzle, and to further explore 
the importance of triaxial deformation in this mass re¬ 
gion, we apply the cranked self-consistent Hartree-Fock- 
Bogoliubov (CHFB) method and TPSM to the rotational 
properties of neutron-rich, even-even Mo and Ru iso¬ 
topes. 

This article is organized as follows. Section [IT] gives 
a brief introduction to CHFB and TPSM models used 
in this work. In Sec. m the results of calculations for 
Mo and Ru isotopes are presented and compared with 
experiment. Therein, we discuss potential energy and 
routhian surfaces, quasi-particle routhian spectra, and 
equilibrium deformations. To test the stability of CHFB 
minima with respect to angular momentum orientation, 
we carry our tilted-axis cranking calculations employing 
the Kerman-Onishi conditions. Finally, the conclusions 
of this work are given in Sec. |IV[ 


II. THE MODEL 

A. Cranked Skyrme-Hartree-Fock-Bogoliubov 
Model 

Our CHFB calculations were performed with the DFT 
solver HFODD (version 2.49t) [33]. Parity and y- 
signature {Ry = exp (—itt Jj,)) symmetries are conserved; 
the corresponding quantum numbers are denoted as tt 
and r. The quasi-particle HFB wave functions were ex¬ 
panded in 800 spherical harmonic oscillator basis states 
with the oscillator frequency of hui =49.2 MeV/A^/^. We 
have tested that such a basis provides a very reasonable 
precision for the observables studied. 

In the particle-hole channel, we employ the global 
Skyrme energy density functional UNEDFO optimized in 
Ref. |35| . In the pairing channel, we take the zero-range 
density-dependent pairing force |36| with the Lipkin- 
Nogami correction for particle-number fluctuations. The 
original pairing strengths are taken as (V)f^,V) 7 ) = 

(-170.374,-199.202) MeV fm^, with a cutoff energy in 
the quasi-particle spectrum of Acut = 60 MeV. In the 
present calculation, the strengths of the pairing force for 
neutrons and protons have been increased by 5% to re¬ 
produce the kinematic moment of inertia of the g.s. band 
(g.b.) of ^°®Mo. As discussed below, the calculated po¬ 
tential energy surfaces are not sensitive to such a small 
variation of pairing strengths. 

In the multi-dimensional potential energy surface cal¬ 
culations, the constraints are imposed on expectation 
values of multipole moments. We use the Augmented 
Lagrangian Method m to perform the constrained it¬ 
erations. The total routhians were computed within the 
principal-axis cranking approach [5]. However, to study 
the stability of the resulting triaxial minima with respect 
to the orientation of the angular momentum vector, we 
applied the Kerman-Onishi conditions implemented as 
in Refs. [35] [33]. Since the Lipkin-Nogami method is not 
strictly variational, the g.b. minimum at nonzero angular 
momentum was obtained by minimizing the constrained 
total routhian surface. This increases precision of calcu¬ 
lations, especially when the minima are soft |40| . 


B. Triaxial Projected Shell Model 

Recently, multi-quasiparticle TPSM approach has 
been developed and it has been shown to provide a con¬ 
sistent description of yrast, 7 (A = 2) and 77 {K = 4) 
bands in transitional nuclei [m 133]. In this method, 
the three dimensional projection technique is employed 
to project out the good angular-momentum states from 
product states built upon quasiparticle (q.p.) configura¬ 
tions of triaxially deformed Nilsson-|-BCS model. The 
shell model Hamiltonian is then diagonalized in this 
angular-momentum projected basis. The TPSM space 
includes multi-quasiparticle states; hence, is capable 
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of describinff near-yrast band structures at hiffh-spins 

m^m- 

The TPSM basis employed in this study consists of 0- 
q.p. vacuum, two-proton, two-neutron, and the four-q.p. 
configurations [iS]. The q.p. basis chosen is adequate 
to describe high-spin states up to angular momentum 
/ ^ 20. In the present analysis we shall, therefore, re¬ 
strict our discussion to this spin regime. As in the earlier 
TPSM calculations, we use the pairing plus quadrupole- 
quadrupole Hamiltonian [ 331117 ]: 

fj, 

where Hq is the single-particle spherical Nilsson Hamil¬ 
tonian, X is the strength of the quadrupole-quadrupole 
force related in a self-consistent way to deformation of 
the q.p. basis, and Gm and Gq are the strengths of 
the monopole and quadrupole pairing terms, respectively. 
The configuration space employed corresponds to three 
principal oscillator shells Afosc- i^[3,4, 5] and 7 r[ 2 , 3,4]. 
The pairing strengths have been parametrized as in 
Refs. |3UI35| in terms of two constants Gi and G 2 - In this 
work, we choose Gi = 16.22 MeV and G 2 = 22.68 MeV; 
with these pairing strengths we approximately reproduce 
the experimental odd-even mass differences in this re¬ 
gion. The quadrupole pairing strength Gq is assumed 
to be proportional to Gm, and the proportionality con¬ 
stant was set to 0.18 |32l|48|. The single-particle basis is 
that of the deformed Nilsson Hamiltonian parametrized 
in terms of axial (e) and triaxial (e') deformations re¬ 
lated to the standard Bohr triaxiality parameter 7 by 
7 = arctan(£'/e). 

III. RESULTS AND DISCUSSIONS 
A. CHFB Results 

1. Ground-state Potential Energy Surfaces 

The g.s. UNEDFO PESs for loe.iosjyjg 
i°®’110412Ru are shown in Eigs.[^and[^ respectively. All 
these nuclei are predicted to be triaxial in their ground 
states. It is seen that the PESs are practically not af¬ 
fected by a 5% change in the pairing strengths. In par¬ 
ticular, the triaxial minima appearing at {Q 20 ,Q 22 ) ~ 
(8.0 — 9.5, 2.0 — 3.0) b are not affected by pairing. The 
corresponding g.s. quadrupole deformations (/? 2 ; 7 ) are 
displayed in Table|l] For predict the triaxial 

g.s. minima at (/32,7) ~ (0.19,17°). Similar results were 
obtained in the macroscopic-microscopic calculations of 
Refs. [T] Uni HH] and HFB-I-DIS calculations [22]. For 
i 08 ,iio,ii 2 j^^^ we also predict triaxial g.s. minima; this 
is consistent with Refs. [Zl|3l] and HF-I-SHI calculations 
of Ref. [T2|. Triaxial g.s. shapes for were also 

obtained in the survey | 1 S| but ^^^Ru was calculated to 
be axial. 



FIG. 1. (Color online) PES in (Q 20 ,Q 22 ) plane in 
CHFB+UNEDFO for i°®Mo and i°®Mo. Left: standard pair¬ 
ing strengths. Right: pairing strengths increased by 5%, see 
text. The difference between contour lines is 0.5 MeV. 



Q20(b) 


FIG. 2. (Color online) Similar to Fig.[^ but for ^°®Ru, ^^°Ru, 
and 


8. High spin behavior 

The angular momentum alignment pattern of Mo and 
Ru nuclei is governed by the vhii /2 and TTgQ /2 high- 
j q.p. excitations, which give rise to strong shape 
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TABLE I. Bohr quadrupole deformation parameters j3 
and 7 calculated in CHFB+UNEDFO for and 



i“®Mo 

i“®Mo 

i“«Ru 

“ORu 

“"Ru 

02 

0.19 

0.18 

0.16 

0.16 

0.15 

7 

16° 

18° 

24° 

25° 

24° 


polarization effects [?]• Figures and show self- 
consistent CHFB-I-UNEDFO 1-q.p- routhian diagrams 
versus rotational frequency for ^°®Mo and ^^^Ru, respec¬ 
tively. In both cases, the alignment of and 

■^( 99 / 2 )^ pairs occurs at similar rotational frequencies of 
hw «0.3MeV. At higher rotational frequencies, a tran¬ 
sition is expected from the g.b. configuration to aligned 
^^(^ 11 / 2 )^ and 7r((7g/2)^ 2-q.p. bands, and then to a 4-q.p. 
^(^ii/ 2 )^’’‘( 59 / 2 )^ band. These two consecutive cross¬ 
ings are difficult to follow in CHFB calculations, as this 
would require a diabatic-configuration extension [50H52] 
of the current framework. Such an extension is highly 
nontrivial in CHFB as the self-consistent mean-fields as¬ 
sociated with aligned configurations are expected to dif¬ 
fer significantly from those of the g.b. [7]. Moreover, 
pairing correlations in the aligned bands are quenched 
and this causes numerical instabilities around the band 
crossing. Therefore, to provide interpretation of the tran¬ 
sition quadrupole moments at higher angular momenta, 
we carry out cranked Skyrme-Hartree-Fock (CHF) cal¬ 
culations without pairing at hio > 0.3 MeV. In this case, 
diabatic configurations can be defined by the number 
of single-particle routhians occupied in the four parity- 
signature blocks [53]. Specifically, each neutron and pro¬ 
ton configurations is defined by four occupation numbers 
[iV++,iV+_,iV_+, N _] representing the number of par¬ 

ticles IV,r,r occupying single-particle states of given tt and 
r. The lowest total routhian with tt = -|- and r = 1 is 
associated with the yrast configuration. 

The angular momentum alignment (total angular mo¬ 
mentum as a function of rotational frequency) is shown in 
Fig. [^for ^°®Mo, ^°®Ru, and ^^^Ru. Below the predicted 
band crossing at fiw ^0.3 MeV, our calculations repro¬ 
duce experiment well. (Note that our pairing strengths 
were adjusted to match the kinematic moment of inertia 
of ^°®Mo.) The first band crossing, associated with the 
alignment of the / 111/2 neutron pair, is seen in ^°®’^^^Ru 
data slightly below huj = 0.4 MeV, and it is significantly 
delayed in ^°®Mo. The predicted aligned configuration 
above the band crossing has a fairly different shape as 
compared to that of the g.b., and it is difficult to follow 
the g.b. at huj > 0.3 MeV. 

To investigate the evolution of nuclear shapes with ro¬ 
tation, we compute the equilibrium j32 and 7 deforma¬ 
tions for low-lying ir = +,r = 1 bands in 
(Fig. 1^ and 108 , 110 ,ii 2 j^^ (Fig. [^. In all cases consid¬ 
ered, the triaxial paired g.b. undergoes small centrifugal 
stretching in the direction of / 32 . For instance, in the case 



0.0 0.1 0.2 0.3 0.0 0.1 0.2 0.3 

hoi (MeV) 


FIG. 3. (Color online) One-quasiparticle routhian diagram 
for ^“®Mo obtained with CHFB-I-UNEDFO. The parity tt and 
signature r of individual levels are indicated in the following 
way: tt = +,r = +i - solid line; tt = +,r = —i - dotted line; 
TT = —,r = +i - dot-dashed line; tt — —,r = —i - dashed line. 
The thin line indicates the Fermi energy. 
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FIG. 4. (Color online) Similar to Fig. but for ^^^Ru. 


of ^°®Ru, /32 increases from the value of 0.15 at huj = Q 
to 0.17 at huj = 0.3 MeV. 

At higher spins (10 < I < 36), when pairing is ne¬ 
glected in our calculations, it is useful to label many- 
body configurations [iV_|_+, _ ,N _V_] by the num¬ 

ber of occupied intruder levels, i.e., Afosc = 4 protons (pri¬ 
marily 59 / 2 ) and Afosc = 5 neutrons (primarily ft-n/g)- 
For instance, the aligned configuration 7r(9,9,12,12) 0 
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FIG. 5. (Color online) Angular momentum alignment for 
“®Mo and CHFB (fio; < 0.3 MeV) and CHF 

{hu) > 0.3 MeV) calculations are compared to experiment 

[I1HI6] 


1/(17,17,15,15) in (shown by circles in Fig. [^a)) 

can be denoted as and the same label applies to 

the 7 r( 9 ,9,12,12) 0 ^(18,18,15,15) in ^°®Mo (shown by 
up-triangles in Fig. |^b)). 

The quadrupole deformations /32 of aligned bands are 
predicted to be in the range of 0.12 < [32 < 0.16, which 
represents a reduction as compared to the shapes of 
paired ground-state bands. The aligned bands remain 
triaxial with 7 values around —30° up to huj = 0.6 MeV. 
This finding is consistent with the results of early Ref. [7] 
employing cranked macroscopic-microscopic approach. 
At the highest rotational frequencies considered, our cal¬ 
culations predict the appearance of aligned triaxial con¬ 
figurations with 7 > 0 , which eventually terminate at 
oblate shapes (7 « 60°), see, e.g.. Fig. [^b). 

To study the stability of different triaxial minima at 
high spins, we analysed related diabatic total routhians in 
the (<520: Q 22 ) plane. In Fig.|^we show the total routhian 
surfaces at huj = 0.5 MeV for the selected low-lying 
aligned configurations in ^°®Mo discussed in Fig. [^a). 
For all those configurations, the collective triaxial mini¬ 
mum with 7 between —30° and —15° appears lowest in 
energy. For the configuration 7 r 4 ^i/ 5 ^ shown in Fig. |^a), 
we also predict a noncollective oblate state with / = 34 
that represents a termination point of 7 > 0 band. 

To eliminate spurious minima that are unstable with 
respect to the angular momentum orientation, we also 
investigated the dependence of the routhians on the an¬ 
gular momentum tilting angle 9 with respect to the axis 
of rotation (y-axis). To this end, we used the tilted- 
axis-cranking formalism of Refs. [351ISH]- The calcula¬ 
tions were performed for the aligned bands in ^°®Mo At 




P cos(y+30°) 


FIG. 6. (Color online) Summary of equilibrium deformations 
of the lowest tt = -|-,r = 1 bands in calculated 

with CHFB-fUNEDFO (ground band) and CHF-fUNEDFO 
(aligned bands). The rotational frequency is varied from zero 
to Two = 0.6 MeV. The corresponding range of angular momen¬ 
tum is indicated. The aligned band are classified according to 
the number of occupied high-A/” intruder levels (A/” = 5 and 4 
for neutrons and protons, respectively). 


hut < 0.5 MeV, the total routhians of triaxial (7 < 0) 
configurations 7 r 4 ®j/ 5 ®, and 7 r 4 ®i/ 5 '^ of Fig. [^a) 

show a minimum at d = 0°. At feu « 0.5 MeV, the 
routhians become very soft in 9, indicating a large- 
amplitude collective motion in this direction. This in¬ 
stability is not present for a (tt = — ,r = 1 ) configu¬ 
ration 7 r( 9 , 9,12,12) G) :/(18,17,15,14) ( 7 r 4 ^i/ 5 ®), which 
shows a pronounced minimum at 0 = 90° associated with 
7 > 0. This result is consistent with the deformation¬ 
driving effect of aligned /iii /2 neutrons orbitals discussed 
in Ref. |7]. 

The transition quadrupole moments along the yrast 
band in lOFios^o and i 08 .iio,ii 2 j^^ 

as a function of rotational frequency. At low rotational 
frequencies fiw < 0.3 MeV, there is a gradual increase of 
Qt with uj due to the centrifugal stretching effect seen 
in Figs. and As discussed earlier, at higher frequen¬ 
cies cranking calculations are performed without pairing. 
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FIG. 7. (Color online) Similar to Fig.[^ but for 

While this approximation seriously affects the predicted 
angular momentum alignment shown in Fig.[^ the equi¬ 
librium shapes obtained in the CHF are reasonable ap¬ 
proximations to those obtained in the full CHFB frame¬ 
work , and reproduce experimental Qt-values for 



FIG. 8. (Color online) Diabatic total routhian surfaces for 
“®Mo at huj = 0.5 MeV calculated in CHF-bUNEDFO for the 
configurations: (a) 7r(9, 9,12,12) G) 1/(17,17,15,15) 
in Fig. [^a)); (b) 7r(10,10,11,11) (g) j/( 17, 17,15,15) 
and (c) 7r(10,10,11,11) (g) z/(16,16,16,16) (7r4®t/5®). 

aligned configurations [57l[58]. As seen in Fig.|^ the pre¬ 
dicted transition quadrupole moments in aligned bands 
are slightly reduced with respect to the low-spin region 
due to the deformation reduction associated with the 
aligned J^/iii /2 and 7159/2 pairs. This reduction is gen¬ 
erally consistent with experiment, except perhaps for 
^^°Ru, where theory overestimates the measured Qt val¬ 
ues above fiw = 0.3 MeV. 


B. TPSM Results 

The TPSM calculations proceed in several stages. In 
the first stage, the deformed basis is constructed based 
on the eigenstates of the triaxially deformed Nilsson po¬ 
tential. The (32 deformation has been chosen such that 
the lowest quadrupole transition from 2 + —> 0 "*" is repro¬ 
duced. The non-axial deformation parameter 7 is chosen 
from the minimum of the g.s. PES obtained in TPSM 
calculations. For ^°®Mo and where the PES 

minima are 7 -soft, triaxial deformation was adjusted to 
the experimental bandhead energy of the 7 -band as it is 
known to be quite sensitive to the 7 -deformation. The 
adopted values of 7 are listed in Table |ll] It is seen that 
strongly triaxial shapes are expected in all cases, and this 
confirms the CHEB-I-UNEDEO results. 

In the next step, the good angular-momentum basis is 
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FIG. 9. (Color online) Transition quadrupole moments for 
i06,108mo and i 08 ,iio,ii 2 j^^ calculated in CHFB (open circles) 
and CHF (dots) compared to experiment. The Qt value at 
7 = 2 is taken from Ref. and the high-spin values from 
Ref. [23]. 


TABLE H. Triaxial quadrupole deformation parameters 
7 employed in the TPSM calculation for and 

108.110,112Ru isotopes. 


Nucleus 


i°®Mo 

i°«Ru 

’i°Ru 

“"Ru 

7 

to 

0 

0 

25° 

29° 

to 

00 

0 

25° 


obtained from the triaxial Nilsson-l-BCS wave functions 
by using the three-dimensional angular-momentum pro¬ 
jection operator. Finally, the shell model Hamiltonian 
Q is diagonalized in this good-angular-momentum ba¬ 
sis. 


1. Band structures 


The TPSM band structures and the associated elec¬ 
tromagnetic transition rates obtained in TPSM are quite 
rich and will be discussed in a separate paper [50] . Figure 
Fig. 10 shows the results for as an representative 

example, as in this nucleus all bands are known up to high 
spins. It is seen that TPSM reproduces the experimental 
band energies quite well. The calculations slightly over- 


“ 2 Ru 


g-band 


y-band 


YY-band 



^ ^ TPSM Expt TPSM Expt TPSM Expt 

FIG. 10. (Color online) Comparison between experimental 
|59 | and calculated band structure of ”^^Ru. 


estimate the bandhead energy of the 77 -band; a similar 
result was also obtained for, e.g., ^°®Mo and ^°®Ru. The 
first q.p. ( 1111 / 2 )^ neutron alignment is predicted around 
1 = 8, and the transition to a 4-q.p. nu(/iii/ 2 )^ 7 r(gg/ 2 )^ 
band is expected to occur around I = 16. 



FIG. 11. (Color online) Comparison of experimental and cal¬ 
culated band structures in ’^'’®Mo at four values of triaxial 
deformation 7 . 

To illustrate the importance of the 7 degree of freedom 
in the description of the band structures in the Mo- and 















































Ru-isotopes, we have carried out TPSM calculations for 
lo^Mo for a range of 7 -values. The obtained band struc¬ 
tures at 7 = 5°, 15°, 20°, and 25° are shown in Fig.[^ At 
7 = 5°, the calculated 7 - and 77 -bands are shifted with 
respect to experiment by more than IMeV. At 7 = 15°, 
7 - and 77 -bands are shifted down in energy, and it is at 
7 = 20 ° that all the three bands are reproduced fairly 
accurately. At still higher value of 7 = 25°, the devi¬ 
ation from experiment grows again. For ^°®Mo, there¬ 
fore, 7 = 20 ° is the optimum triaxial deformation in 
TPSM. Interestingly, the same value of 7 is predicted 
by CHFB-^UNEDFO, see Table 

2. Transition Quadrupole Moments 

Using the TPSM wave functions and standard E2 ef¬ 
fective charges (e„ = 0.5e and Cp = 1.5e) we have evalu¬ 
ated the transition quadrupole moments along the yrast 
line of the studied Mo- and Ru-isotopes, see Fig.[^ The 
overall behavior of the measured Qt-values is reproduced 
quite well by the TPSM approach. The drop in Qt ob¬ 
served for all the studied isotopes around / = 8 is due to 
the quasiparticle alignment of the /iii /2 neutron pair, and 
in some isotopes there is a further drop around / = 16 
due to a consecutive alignment involving the 59/2 pro¬ 
ton pair. Figure [T^ a) also shows the TPSM predictions 
for Qt in ^°®Mo at 7 = 5°, 15°, and 25°. Again, it is 
seen that the best reproduction of experimental data is 
obtained at the optimum value of 7 = 20 °. 

IV. CONCLUSIONS 

Stimulated by the recent experimental data on the 
transition quadrupole moments in the rotational bands 
of neutron-rich Mo and Ru nuclei [23] , we studied 
the shapes of the band structures in ^°®’^°®Mo, and 
i 08 ,iio,ii 2 j^^^ We employed two complementary theoret¬ 
ical models: self-consistent CHFB-I-UNEDFO approach 
and TPSM. 

The triaxial PESs obtained in CBFB-I-UNEDFO show 
stable g.s. triaxial minima in all cases. At higher angu¬ 
lar momenta, the consecutive band crossings along the 
yrast line are expected, associated with the alignment 
of s-iid ^( 39 / 2 )^ pairs. The quadrupole defor¬ 

mations (^2 in the aligned bands are predicted to be re¬ 
duced, but the shapes remain strongly triaxial. This 
result confirms predictions of an earlier works [71 in] 
based on the cranked macroscopic-microscopic method 
with the Woods-Saxon average potential. The decrease of 
the corresponding transition quadrupole moments above 
fujj ^ 0.3 MeV reflects the change in ^2 and 7 due to 
the q.p. alignment, and this reduction is consistent with 
experiment. 

The results obtained with TPSM paint the same pic¬ 
ture as CBFB calculations, and strongly favor the triaxial 
interpretation. Indeed, both the energies of yrast, 7 , and 



FIG. 12. (Color online) Transition quadrupole moments for 
i06,108mo and i 08 .iio.ii 2 j^^ calculated in TPSM (solid line) 
compared to experiment. The Qt value at / = 2 is taken 
from Ref. |53| and the high-spin values from Ref. |23| . For 
^°®Mo, we also show TPSM results at 7 = 5°, 15°, and 25° 
(dash-dotted lines). 


77 bands, and transition quadrupole moments are well 
described assuming stable triaxial shapes. Similar as in 
our CBFB calculations, transition quadrupole moments 
obtained in TPSM exhibit a reduction at neutron and 
proton band crossings. 

In summary, according to our analysis, high-spin be¬ 
havior of ^°®’^°®Mo, and i 08 .iio.ii 2 j|^ jg consistent with 
triaxial rotation. The predicted triaxial g.s. minima are 
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fairly shallow, and this perhaps is why in some calcu¬ 
lations, e.g., the cranked relativistic Hartree-Bogoliubov 
model of Ref. [53], axial configurations may be slightly 
favored. 


ACKNOWLEDGMENTS 

Useful discussions with Jacek Dobaczewski and Nico¬ 
las Schunck are gratefully acknowledged. This material 
is based upon work supported by the U.S. Department of 
Energy, Office of Science, Office of Nuclear Physics un¬ 
der Award Numbers No. DOE-DE-SC0013365 (Michigan 
State University) and No. DE-SC0008511 (NUCLEI Sci- 
DAC Collaboration). An award of computer time was 
provided by the National Institute for Computational 
Sciences (NICS) and the Innovative and Novel Compu¬ 
tational Impact on Theory and Experiment (INCITE) 
program using resources of the OLCF facility. 


[1] J. Eberth, R. A. Meyer, and K.Sistemich, eds.. Nuclear 
Structure of the Zirconium Region (Springer, 1988). 

[2] J. Stachel, N. Kaffrell, E. Grosse, H. Emling, H. Folger, 

R. Kulessa, and D. Schwalm, Nucl. Phys. A 383, 429 
(1982) 

[3] J. Srebrny, T. Czosnyka, and C. Droste, Nucl. Phys. A 
766,25 (2006) 

[4] Y. Luo, S. Zhu, J. Hamilton, J. Rasmussen, A. Ramayya, 
C. Goodin, K. Li, J. Hwang, D. Almehed, S. Frauendorf, 

V. Dimitrov, J. ye Zhang, X. Ghe, Z. Jang, I. Stefanescu, 
A. Gelberg, G. Ter-Akopian, A. Daniel, M. Stoyer, R. Do- 
nangelo, J. Gole, and N. Stone, Phys. Lett. B 691, 285 
( 2010 ) 

[5] Wrzosek-Lipska, Phys. Rev. C 86, 064305 (2012) 

[6] A. G. Smith, J. L. Durell, W. R. Phillips, W. Urban, 
P. Sarriguren, and I. Ahmad, Phys. Rev. G 86, 014321 
( 2012 ) 

[7] J. Skalski, S. Mizutori, and W. Nazarewicz, Nucl. Phys. 
A 617, 282 (1997) 

[8] A. Bohr and B. R. Mottelson, Nuclear Structure, vol. II 
(W. A. Benjamin, Reading, 1975). 

[9] S. Frauendorf, Rev. Mod. Phys. 73, 463 (2001) 

[10] E. Cheifetz, R. G. Jared, S. G. Thompson, and J. B. 
Wilhelmy, Phys. Rev. Lett. 25, 38 (1970) 

[11] M. Hotchkins, J. Durell, J. Fitzgerald, A. Mowbray, 

W. Phillips, I. Ahmad, M. Carpenter, R. Janssens, 
T. Khoo, E. Moore, L. Morss, P. Benet, and D. Ye, 
Nucl. Phys. A 530, 111 (1991) 

[12] Aysto, P. Jauho, Z. Janas, A. Jokinen, J. Parmonen, 
H. Penttila, P. Taskinen, R. Beraud, R. Duffait, A. Em- 
sallem, J. Meyer, M. Meyer, N. Redon, M. Leino, K. Es- 
kola, and P. Dendooven, Nucl. Phys. A 515, 365 (1990) 

[13] J. Shannon, W. Phillips, J. Durell, B. Varley, W. Urban, 
C. Pearson, I. Ahmad, C. Lister, L. Morss, K. Nash, 
C. Williams, N. Schulz, E. Lubkiewicz, and M. Bentaleb, 
Phys. Lett. B 336, 136 (1994) 

[14] Q. H. Lu, K. Butler-Moore, S. J. Zhu, J. H. Hamilton, 
A. V. Ramayya, V. E. Oberacker, W. C. Ma, B. R. S. 
Babu, J. K. Deng, J. Kormicki, J. D. Cole, R. Aryaeine- 
jad, Y. X. Dardenne, M. Drigert, L. K. Peker, J. O. Ras¬ 
mussen, M. A. Stoyer, S. Y. Chu, K. E. Gregorich, 1. Y. 
Lee, M. F. Mohar, J. M. Nitschke, N. R. Johnson, F. K. 
McGowan, G. M. Ter-Akopian, Y. T. Oganessian, and 

J. B. Gupta, Phys. Rev. C 52, 1348 (1995) 

[15] H. Hua, C. Y. Wu, D. Cline, A. B. Hayes, R. Teng, R. M. 


Clark, P. Fallon, A. Goergen, A. O. Macchiavelli, and 

K. Vetter, Phys. Rev. C 69, 014317 (2004) 

[16] C. Y. Wu, H. Hua, D. Cline, A. B. Hayes, R. Teng, D. Ri¬ 
ley, R. M. Clark, P. Fallon, A. Goergen, A. O. Macchi¬ 
avelli, and K. Vetter, Phys. Rev. C 73, 034312 (2006) 

[17] J. H. Hamilton, S. J. Zhu, Y. X. Luo, A. V. Ramayya, 

J. O. Rasmussen, J. K. Hwang, J. G. Wang, Q. Xu, H. B. 
Ding, S. H. Liu, G. M. Ter-Akopian, A. V. Daniel, and 
Y. Oganessian, AlP Conf. Proc. 1175, 166 (2009) 

[18] J. H. Hamilton, S. J. Zhu, Y. X. Luo, A. V. Ramayya, 

S. Frauendorf, J. O. Rasmussen, J. K. Hwang, S. H. Liu, 
G. M. Ter-Akopian, A. V. Daniel, and Y. Oganessian, 
Nucl. Phys. A 834, 28c (2010) 

[19] A. G. Smith, J. L. Durell, W. R. Phillips, M. A. 
Jones, M. Leddy, W. Urban, B. J. Varley, 1. Ahmad, 

L. R. Morss, M. Bentaleb, A. Guessous, E. Lubkiewicz, 
N. Schulz, and R. Wyss, Phys. Rev. Lett. 77, 1711 
(1996) 

[20] W. Urban, T. Rzaca-Urban, J. Durell, W. Phillips, 
A. Smith, B. Varley, I. Ahmad, and N. Schulz, Eur. 
Phys. J. A 20, 381 (2004) 

]21] H. Watanabe, K. Yamaguchi, A. Odahara, T. Sumikama, 
S. Nishimura, K. Yoshinaga, Z. Li, Y. Miyashita, and 

K. Sato, Phys. Lett. B 704, 270 (2011) 

[22] S. Zhu, Y. Luo, J. Hamilton, J. Rasmussen, A. Ramayya, 

J. Hwang, H. Ding, X. Che, Z. Jiang, P. Gore, E. Jones, 

K. Li, I. Lee, W. Ma, G. Ter-Akopian, A. Daniel, 

S. Frauendorf, V. Dimitrov, J. Zhang, A. Gelberg, I. Ste- 
fancescu, and J. Cole, Prog. Part. Nucl. Phys. 59, 329 
(2007) 

[23] J. Snyder, W. Reviol, D. Sarantites, A. Afanasjev, 
R. Janssens, H. Abusara, M. Carpenter, X. Chen, 

C. Chiara, J. Greene, T. Lauritsen, E. McCutchan, 

D. Seweryniak, and S. Zhu, Phys. Lett. B 723, 61 
(2013) 

[24] P. Moller, R. Bengtsson, B. G. Carlsson, P. Olivius, and 

T. Ichikawa, Phys. Rev. Lett. 97, 162502 (2006) 

[25] Moller, R. Bengtsson, B. Carlsson, P. Olivius, 
T. Ichikawa, H. Sagawa, and A. Iwamoto, At. Data Nucl. 
Data Tab. 94, 758 (2008) 

[26] M. Boyiikata, P. V. Isacker, and I. Uluer, J. Phys. G 37, 
105102 (2010) 

[27] K. Nomura, N. Shimizu, and T. Otsuka, Phys. Rev. C 
81, 044307 (2010) 

[28] I. Bentley and S. Frauendorf, Phys. Rev. C 83, 064322 




10 


( 2011 ) 

[29] R. Rodriguez-Guzman, P. Sarriguren, L. Robledo, and 
S. Perez-Martin, Phys. Lett. B 691, 202 (2010) 

[30] F. R. Xu, P. M. Walker, and R. Wyss, Phys. Rev. C 65, 
021303 (2002) 

[31] J. Q. Faisal, H. Hua, X. Q. Li, Y. Shi, F. R. Xu, H. L. 
Liu, Y. L. Ye, and D. X. Jiang, Phys. Rev. C 82, 014321 
( 2010 ) 

[32] A. Bhat, A. Bharti, and S. K. Khosa, Int. J. Mod. Phys. 
E 21, 1250030 (2012) 

[33] Y.-X. Liu and Y. Sun, J. Phys. Conf. Ser. 420, 012046 
(2013) 

[34] N. Schunck, J. Dobaczewski, J. McDonnell, 
W. Nazarewicz, J. Sheikh, A. Staszczak, M. Stoitsov, 
and P. Toivanen, Comput. Phys. Commun. 183, 166 
( 2012 ) 

[35] M. Kortelainen, T. Lesinski, J. More, W. Nazarewicz, 
J. Sarich, N. Schunck, M. V. Stoitsov, and S. Wild, 
Phys. Rev. C 82, 024313 (2010) 

[36] J. Dobaczewski, W. Nazarewicz, and M. V. Stoitsov, 
The Nuclear Many-Body Problem 2001, edited by 
W. Nazarewicz and D. Vretenar (Springer, 2002) p. 181. 

[37] A. Staszczak, M. Stoitsov, A. Baran, and 
W. Nazarewicz, Eur. Phys. J. A 46, 85 (2010) 

[38] Y. Shi, J. Dobaczewski, S. Frauendorf, W. Nazarewicz, 
J. C. Pei, F. R. Xu, and N. Nikolov, Phys. Rev. Lett. 
108, 092501 (2012) 

[39] Y. Shi, C. L. Zhang, J. Dobaczewski, and 
W. Nazarewicz, Phys. Rev. C 88, 034311 (2013) 

[40] M. V. Stoitsov, J. Dobaczewski, W. Nazarewicz, S. Pittel, 
and D. J. Dean, Phys. Rev. C 68, 054312 (2003) 

[41] J. A. Sheikh and K. Hara, Phys. Rev. Lett. 82, 3968 
(1999) 

[42] Y. Sun, K. Hara, J. A. Sheikh, J. G. Hirsch, V. Velazquez, 
and M. Guidry, Phys. Rev. C 61, 064323 (2000) 

[43] J. A. Sheikh, Y. Sun, and R. Palit, Phys. Lett. B 507, 
115 (2001) 

[44] Y. Sun, J. A. Sheikh, and G.-L. Long, Phys. Lett. B 
533,253 (2002) 


[45] J. A. Sheikh, G. H. Bhat, Y. Sun, G. B. Vakil, and 
R. Palit, Phys. Rev. C 77, 034313 (2008) 

[46] J. Sheikh, G. Bhat, R. Palit, Z. Naik, and Y. Sun, Nucl. 
Phys. A 824, 58 (2009) 

]47] K. Hara and Y. Sun, Int. J. Mod. Phys. E 04, 637 (1995) 

[48] E. Y. Yeoh, S. J. Zhu, J. H. Hamilton, K. Li, A. V. Ra- 
mayya, Y. X. Liu, J. K. Hwang, S. H. Liu, J. G. Wang, 

Y. Sun, J. A. Sheikh, G. H. Bhat, Y. X. Luo, J. O. Ras¬ 
mussen, 1. Y. Lee, H. B. Ding, L. Gu, Q. Xu, Z. G. Xiao, 
and W. C. Ma, Phys. Rev. C 83, 054317 (2011) 

[49] D. Galeriu, D. Bucurescu, and M. Ivascu, Journal of 
Physics G: Nuclear Physics 12, 329 (1986) 

[50] S. Cwiok, W. Nazarewicz, J. Dudek, J. Skalski, and 

Z. Szymahski, Nucl. Phys. A 333, 139 (1980) 

[51] Z. Szymariski, Fast nuclear rotation (Oxford University 
Press, 1983). 

[52] A. Axelsson, R. Bengtsson, and J. Nyberg, Nucl. Phys. 
A 708, 245 (2002) 

[53] J. Dobaczewski and J. Dudek, Comput. Phys. Commun. 
131, 164 (2000) 

]54] S. Raman, C. Malarkey, W. Milner, C. N. Jr., and 
P. Stelson, Atomic Data and Nuclear Data Tables 36, 
1 (1987) 

[55] W. Nazarewicz, J. Dudek, R. Bengtsson, T. Bengtsson, 
and 1. Ragnarsson, Nucl. Phys. A 435, 397 (1985) 

[56] J. Dudek, W. Nazarewicz, and N. Rowley, Phys. Rev. C 
35, 1489 (1987) 

]57] W. Satula, J. Dobaczewski, J. Dudek, and 
W. Nazarewicz, Phys. Rev. Lett. 77, 5182 (1996) 

[58] M. Matev, A. V. Afanasjev, J. Dobaczewski, G. A. 
Lalazissis, and W. Nazarewicz, Phys. Rev. G 76, 034304 
(2007) 

]59] C. Xing-Lai, Z. Sheng-Jiang, J. H. Hamilton, A. V. 
Ramayya, J. K. Hwang, J. O. Rasmussen, Y. X. Luo, 
C. Yong-Jing, L. Ming-Liang, D. Huai-Bo, U. Yong-Nam, 
1. Y. Lee, and W. C. Ma, Ghinese Phys. Lett. 23, 328 
(2006) 

[60] G. H. Bhat and J. A. Sheikh, to be published (2015). 







